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SUMMARY 


Studies of rotor and propeller performance and induced potential 
flcwfields are made on the basis of a rotating actuator disk concept, 
with special emphasis placed on rotors hovering out of ground effect. 

A new theory for the optimum performance of rotors hovering 0(22 is 
developed and presented. An extended theory for the optimum performance 
of rotors and propellers in axial motion is also presented. Numerical 
results are presented for the optimum distributions of blade-bound cir- 
culation together with axial inflow and ultimate wake velocities for 
the hovering rotor over the range of thrust coefficient of interest in 
rotorcraft applications. Shapes of the stream tubes and of the velocities 
in the slipstream are obtained, using available methods, for optimum and 
"off -optimum" circulation distributions for rotors hovering in and out of 
ground effect. Results of the optimum performance study show that 
improvements made possible by the present performance theory, which fully 
account for the effects of slipstream rotation, is important to the 
optimum circulation as well as to the power and the thrust requirements 
of the rotor in the higher thrust coefficient range. A number of explicit 
formulae useful in computing rotor and propeller induced flows are pre- 
sented for stream functions and velocities due to distributions of cir- 
cular vortices over axi-symmetric surfaces. The computed slipstream 
shape shows that the assumption of a uniform contraction ratio is 
reasonable and that, within the limitations of the infinitely bladed model, 
the computed results provide a reasonable description of the wake boundary. 
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NOMENCLATURE 


a 

c 


T 


D 

-♦ 

e 


E 

F 

— • 

F 

G 

G 

] 

h 

H 

k 



“2 

M 


radial location 01 ' perturbation in the ultimate wake 
contraction ratio 

power coefficient defined in Eq. (56) 

thrust coefficient defined in Eq. (55) 

domain of integration consisting of the rotor slipstream 

unit vector 

complete elliptic integral of the second kind 
function defined in Eq. (B-10) 
body force vector 

Green's function for the streamfunction defined by Eq. (23) 

Green's function for the axial velocity defined by Eq. (77) 

function defined by Eqs. (43) and (B-6) 

total head of the fluid 

modulus of elliptic integrals 

modulus defined by Eq. (A-26) 

modulus defined by Eq. (A -27 ) 

modulus defined by Eq. (A-54) 

modulus defined by Eq. (A-55) 

complete elliptic integral of the first kind 

indicates the segment of constant circulation 

total number of segments of constant circulation 

characteristic defined by Eq. (A-51) 

characteristic defined by Eq. (A-52) 

figure of merit defined by Eq. (54) 
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n the "characteristic" of elliptic integrals of the third kind 

n unit outward normal vector to the control volume 

N a constant defined by Eq. (45 ) 

p static pressure of the fluid 

P power expended by the rotor 

q velocity vector 

r radial coordinate in a cylindrical coordinate system 

r position vector of an observation point 

r^ inner radius of an annular surface 

r 2 outer radius of an annular surface 

r' position vector of a vortex filament 

R radial location of a vortex filament 

radius of the rotor or propeller disk 

R^ radial location of the l discontinuity in disk circulation 

B w radius of the slipstream in the ultimate wake 

ft radius of a control volume centered at the rotor hub 

s annular element of area 

s arc length of a vortex segment 

the portion of the control surface which is cut out by the 
wake slipstream 

S 2 the portion of the control surface excluding 

t time 

T thrust of the rotor 

u radial component of the velocity 

v tangential component of the velocity 

w axial component of the velocity 
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W constant value of the axial velocity in the ultimate wake within 

a streamtube of constant circulation 

2 

x is defined as 

2 axial coordinate in a cylindrical coordinate system 

Z axial location of a vortex filament 

a variable of integration of elliptic integrals 

y vortex strength per unit length 

T bound circulation at the disk 

6 function defined by Eq. (A- 12) 

A indicates a perturbed value in the ultimate wake 

e perturbed value of the angular velocity in the ultimate wake 

£ axial component of vorticity 

T) tangential component of vorticity 

9 tangential coordinate in a cylindrical coordinate system 

x vor ex filament strength 

X vortex strength per unit length 

X I constant in Eq. (A -43 ) defining a linear vortex strength 

Xg constant in Eq. (A -43) defining a linear vortex strength 

Xg constant for vortex strength on an annular surface 

A advance ratio defined as w 

§ radial component of the vorticity 

II complete elliptic integral of the third kind 

p density of the fluid 

Y stream function 

ip vorticity vector 

0 angular velocity of the rotor 
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Superscripts 


' indicates a variable of integration 

~ denotes a dimensionless variable defined by Eq. (49) 

Subscripts 

d denotes a variable at the rotor disk 

f refers to a quantity induced by a circular vortex filament 

m indicates a value on the vortex tube 

w refers to a value in the ultimate wake 

o indicates a variable at the rim of the slipstream in the 

ultimate wake 

® denotes a value in the undisturbed freestream 
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1. INTRODUCTION 


The development of a suitable theory for predicting the flow- 
field induced by a rotor has long been a central problem of rotorcraft 
aerodynamics. Since the flowfield induced by the rotor is ultimately 
responsible for the aerodynamic and dynamic behaviors of the rotorcraft, 
improvements in methods for predicting the flowfield contribute to the 
advancement of rotorcraft design in various ways. 

The case of a rotor in hover is of special importance in a number 
of practical problems. These problems include the familiar operational 
problems, the problem of down wash impingement which leads to a grcund 
boundary layer and particle entrainment, the problem of tail rotor 
behavior when immersed in the slipstream of a main rotor, and the 
problem of rotor performance. The general subject of a rotor in for- 
ward flight Involves, of course, many additional important problems and 
is more complicated than the havering case. It is generally expected 
that a thorough understanding of the latter will precede and then con- 
tribute substantially to research of the former. The hovering case, 
however, is in itself of such complexity that a theory sufficiently 
accurate and yet simple enough to be used as a design tool is yet 
unavailable. The development of such a theory is emphasized in the 
current literature. 

A hovering rotor is a special case of a propeller, with zero 
advance velocity. Marine and aircraft engineers have long beer con- 
cerned with propeller aerodynamics. In the majority of situations of 
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interest to the propeller aerodynamic 1st, however, the propeller-induced 
velocity is small in comparison with the propeller's velocity of advance. 
As a consequence, the propeller aerodynamicist is accustomed to simplifi- 
cations that are not necessarily valid for the hovering rotor problem. 
Propeller theories therefore do not usually go over to rotor theories 
in a straight forward manner and predictions based on generalized 
propeller theories do not always correlate well with experimental data 
for rotors (Ref. 1). 

Of the several problems mentioned earlier, the problem of pre- 
dicting the rotor performance in hover out of ground effect (OGE) is 
particularly important. Rotorcraft are typically required to be able 
to hover OGE when occasion demands . Since the power required for a 
rotor to develop a given amount of thrust is the greatest during hover 
OGE, the payload capability of the rotorcraft is generally limited by 
the hovering performance OGE . Mode: :* rotorcraft characteristically 
have small ratios of payload to gross weight at take off. Consequently, 
erro.s in the predicted hovering performance are amplified in the cal- 
culation of payload capability. For example, with a payload to gross 
weight ratio of l/5 at take off, a 5$ error in the prediction of hover- 
ing thrust means a 25$ error in the predicted payload capability. 

In recent years, a sizable effort has been in progress to 
establish suitably simple, yet sufficiently accurate theories and 
methods for predicting hovering rotor performance and induced flowfields. 
A review of the state of art of the hovering rotor problem is presented 
in Section II of this report. This review shows that, although recent 
progress in this area has been extensive, there still exists a need to 
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improve upon the available theories and methods for modern rotors with 
high rotational speed, solidity, and disk loading. 

The purpose of this report is to present recent research •e.outs 
that serve to supplement the available theories and methods. In parti- 
cular, an improved theory for the optimum performance of hovering rotors 
is presented. This theory is based on the actuator disk concert and 
therefore subject to some of the limitations of the well known general 
momentum theory. The correction factor for finite number of blades 
must be estimated separately, for example, by an extension of Lerbs' 
method (Ref. 2) for heavily loaded propellers. The improved theory, 
however, is more complete in that it fully accounts for the effects of 
slipstream rotation. This improvement is not important in the case of 
a lightly loaded rotor, but has finite contributions to thrust and to 
poorer requirements of a heavily loaded rotor. The improved theory per- 
mits the computation, in a straight forward manner, of the optimum 
distributions of circulation and inflow velocity over the rotor disk. 
These distributions, along with the figure of merit and the optimum 
power coefficient, are presented for several values of the thrust coeffi- 
cient. 

In addition to the optimum performance theory, flow fields induced 
by the infinitely bladed rotor hovering in, as well as out of, ground 
effect are presented for several cases of optimum and off -optimum cir- 
culations. The computation of the flowfields is based on the integral 
formulation of Ref. 3* The numerical procedures are those of Ref. 5. 
Several useful analytical expressions developed during the course of 
this research in connection with the computational approach sure pre- 
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sented in Appendix A of this report. A generalized optimum performance 
theory applicable to rotors (or propellers ) in axial motion is presented 
in Appendix B. 
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II . BACKGROUND 


In recent years, the search for better methods of predicting 
rotor induced flowfield and associated performance in hover yielded 
a copious volume of literature on the subject. Several recent articles 
contain comprehensive bibliographies (Refs. 6 and 7) of the earlier work 
as well as review and reassessment of "classical" theories. The brief 
discussion given below provides a summary account of previous theories 
and present efforts. Only a few pertinent articles are referred to 
here. In many cases, a large number of articles exist which are based 
cn essentially the same approach. Only one representative article is 
then mentioned. 

As mentioned earlier, a hovering rotor is a special case of a 
propeller, with zero advance velocity. In the following discussion, 
however, the term "propeller" is used to describe only the case where 
the velocity 0 ? advance is of primary importance. The term "rotor" is 
used to describe che static propeller or the hovering rotor. 

The majority of recent articles on rotor theory are concerned 
with the important problem of optimum performance. The central task 
in the performance study is the determination of a distribution of the 
rotor inflow velocity, over a given rotor disk, that leads to a minimum 
amount of power expenditure for a given amount of thrust. This infor- 
mation, when used in conjunction with the well known blade -element 
theory, permits the prediction of the performance limit and of the 
required radial distribution of circulation in order to approach this 
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limit. It thus provides a rational basis for blade design. The pre- 
diction of the optimum rotor performance, however, does not necessarily 
require a knowledge of the flowfield away from the rotor disk. In fact, 
the major portion of existing literature on rotor theory avoids the 
evaluating of the entire flowfield and treats the subject of optimum 
inflow velocity and optimum performance by themselves. 

The earliest performance theory was the axial momentum theory 
based on the actuator disk (infinitely-bladed propeller) concept for 
propellers (Refs. 6 and 9). Within the context of this theory, the 
propeller operates without any frictional drag on the blade and induces 
no rotation in the slipstream. Since the flow upstream of the pro- 
peller disk is also irrotational, the tangential velocity component at 
the disk is zero. It follows from the Kutta-Joukowski theorem that 
the disk cannot be subjected to a lift force. Disregarding this incon- 
sistency, two well known conclusions were obtained: first, an optimum 

propr" 1 ler is one with a uniform axial velocity at the propeller disk, 
and second, when applied to a hovering rotor, the minimum pov\ ~ coeffi- 
cient is equal to the thrust coefficient to the 3/2 power divided by 
the squrr^ root of 2. The theory provides no information on the radial 
distribution of circulation over the disk and consequently no indication 
as to how the blades should be designed for optimum performance. 

In general, the energy loss due to the rotational motion in the 
slipstream of a propeller is small. The neglect of the slipstream 
rotation wa3 therefore thought to be justified for the purpose of pre- 
dicting the performance of propellers. The actual performance of pro- 
pellers, however, generally does not meet the expectation of the axial 
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momentum theory even when the various loss factors, not including the 
slipstream rotation, are taken into account. This deficiency is 
explained by the general momentum theory, which incorporates a pro- 
cedure to account for some of the effects of slipstream rotation. One 
of the major conclusions of the general momentum theory is that the 
slipstream rotation, though small in terms of the energy content, exerts 
an important influence on the optimum distribution of inflow velocity 
over the propeller disk. A detailed presentation of the general 
momentum theory is given in Ref. 10 which also contains a presentation 
of efforts directed toward the study of helicopter rotors. The general 
momentum theory does provide information on the optimum distribution 
of circulation over the disk. The solution as given in Ref. 10, how- 
ever, is approximate and only partially accounts for the effects of 
slipstream rotation. Some of the neglected effects are of importance 
in the study of hovering rotors. For example, the conservation of 
angular momentum in the slipstream requires an exchange of energy 
between the rotational mode and the axial mode to accompany the con- 
traction of the slipstream. The effect of this exchange of energy is 
neglected in the theory. For propellers where the velocity of advance 
is large compared to the induced velocity, the contraction is small 
and the neglect of the energy exchange is justified as long as the 
rotational energy is not excessively large. For hovering rotors, the 
ultimate wake cross-section is, to the first order, one half the size 
of the rotor disk. The exchange of energy therefore has a significant 
influence on the optimum distribution of circulation over the propeller 
disk for the situation where the rotational energy is not extremely small. 



For rotor applications, it is therefore desirable to improve on the 
general momentum theory by fully sccounting for the effects of slip- 
stream rotation. 

Until very recently, the theory which represents the current 
state of the art for routine calculation of propeller performance was 
the well-established vortex theory. The basic vortex theory represents 
the wake trailing the blades by a distribution of concentric cylindrical 
vortex sheets. These vortex sheets describe the radial variation of 
the blade circulation and thus imply an infinite number of blades. The 
effect of a finite number of blades is obtained approximately by Betz 
and Prandtl (Ref. 11 ) using a tip loss factor. An improved analysis 
by Goldstein (Ref. 12) represents the trailing vortex sheets in the 
ultimate wake, for optimum performance, by dis rete helicoidal surfaces 
of constant helix angle and moving as rigid surfaces. Goldstein's work 
formed the basis of much of the subsequent analyses by Lock, Theodorsen, 
Lerbs, etc., (Refs. 13, l4, 15, 2 and 17) that proved to be sufficiently 
accurate for predicting propeller performance . The use of the vortex 
theory for predicting rotor performance, however, yielded overly opti- 
mistic results at high disk loadings (Ref. l). 

The continuing requirements of higher forward flight speed and 
larger rotorcraft, coupled with the need to keep the rotor small from 
weight and operational flexibility considerations, led to the develop- 
ment of rotors with higher rotational speed and disk loading. This 
resulted in increased inaccuracy of the vortex theory in predicting the 
rotor performance. Reviews of the vortex theory suggest that the wake 
contraction, the slipstream rotation, and the associated non-uniform 
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inflow are the factors that contribute most to the inability of the 
vortex theory to predict the hovering rotor performance accurately. 

This contention is supported by recent results obtained from the 
numerical computation of the entire rotor-induced flowfield. 

In the numerical methods, typically, each rotor blade is repre- 
sented by a lifting line with piecewise uniform circulation. The 
vorticity in the slipstream is represented by a number of discrete 
vortex filaments whose strengths are related to the discontinuities 
in circulation at the blades . Each vortex filaments is in turn approxi- 
mately by a number of straight segments. If the geometry of the trail- 
ing vortex system is known, then the use of the Biot-Savart Law and the 
blade element theory yields a set of simultaneous equations. The 
solution of this set of equations gives the distribution of blade- 
bound vortex strength and hence also the trailing vortex filament 
strengths. The computations of velocities induced by the complete 
vortex system at all points in space, including those on the rotor disk, 
and of the thrust, the power, and the performance of the rotor are then 
straight forward. 

For some rotor configurations and thrust levels, flow visuali- 
zation studies have yielded useful wake geomecry data. Methods that 
utilize experimentally obtained wake geometries in the numerical Biot- 
Savart Law approach described above are often referred to as the "pres- 
cribed-wake" analyses (Refs. 18, 19 ). For cases where experimental 
data permit the establishment of an accurate model for the wake geometry, 
the prescribed-wake analyses have yielded results in good agreement with 
experiments for overall performance (Ref. 7). For a general application 

9 


r 


of the prescribed wake approach, however, accurate experimented, data 
about the weike geometry must be available for a wide range of rotor 
configurations and thrust levels. 

To remove the empiricism inherent in the prescribed wake approach, 
several researchers have developed methods which generate the wake geo- 
metry as a part of the computation procedure. This approach, sometimes 
referred to as the "free-wake" analysis (Ref. 20), employs an iterative 
or a time -step procedure to find, for a given rotor configuration, an 
arrangement of the vortex system which exists in "force -free equili- 
brium". The free-wake analysis generally incorporates a procedure to 
account for the convection of the vortex segments with the fluid but 
neglects the process of viscous diffusion. The analysis requires a 
relatively large amount of computer time for each combination of rotor 
configuration and thrust level. In order to determine rotor configur- 
ations that would give the optimum perf ormance , a parametric study 
involving a large number of geometries must be studied. The amount of 
computer time required by the free-wake analyses in optimum performance 
studies is therefore extremely large. Also important is the fact that 
’ rially the free-wake analyses do not accurately predict the geometry 
of the wake in the region very near the rotor disk (Ref. 21) where the 
vortex system contributes greatly to the inflow velocity at the disk. 

For these reasons, the present utility of the free-wake analyses for 
deoign purposes is limited. 

One advantage of the numerical Biot-Savart Law approach is that, 
in contrast to the previous theories which answer only the questions of 
perfo-nnance , this approach provides information on the flowfield 



surrounding the rotor as well. It must be recognized that the approach 
in use today is based on the inviscid flow analyses and its ability to 
accurately predict the rotor induced flowfield is subject to important 
limitations. In this regard, it is noted that a new numerical method 
which utilizes a generalized Biot-Savart Law and accounts for the 
process of viscous diffusion of vortices has been developed recently 
(Ref. 22). This method is applicable to the rotor problem. 

In the inviscid numerical approach, each rotor blade is typically 
represented by a lifting line with piecewise uniform circulation. Thus, 
in addition to neglecting the viscous diffusion which tends to spread 
the vortex sheets into volume distributions of vorticity, the inviscid 
approach further lumps the surface distribution of vortices into vortex 
filaments for computational conveniences. A method formulated by T. Y. 
Wu (Ref. 3) for heavily loaded propellers, in contrast, utilizes the 
actuator disk concept and distributes the blade-bound vortex azimuthally 
over the propeller disk in an axi-syrametric manner. The trailing vortex 
sheets are now represented by a volume distribution of vortices in an 
axi -symmetric flow. The kinematics of the flow is formulated in terms 
of Stokes' stream function and the rotational velocity component. By 
the use of a Green's function, a non-linear integral equation for the 
stream function is established, allowing the solution for the stream 
function by successive approximations. The flowfield is thus formulated 
as a boundary value problem with a prescribed arbitrary radial distri- 
bution of circulation at the disk. The integral formulation of T. Y. 

Wu for the calculation of the stream function requires numerical quad- 
rature of a triple integral containing a product of Bessel's functions 


11 



as the integrand. Two of the three ranges of integration are infinite. 
For these reasons the computational effort required, though less than 
the Biot-Savart Law approach, is still considerable. 

Greenberg et al, (Refs. 4 and 5) utilized the integral formulation 
to compute flowfields induced by rotors with prescribed uniform and 
piecewise uniform circulation distributions. The Green's function in 
the integral formulation is recognized in their works an the stream 
function due to a circular vortex filament. By expressing the Green's 
function in terms of a Legendre function, the flowfield computations 
were accomplished reasonably rapidly. Similarly, Chaplin (Ref. 23) and 
Cox (Ref. 24) used the integral formulation in studies of axi-symmetric 
flows. 

The actuator disk concept enjoys a substantial advantage in its 
relative simplicity. Results obtained on the basis of the actuator 
disk concept have been successfully used in the past as basic building 
blocks for the understanding of the practical problem of a rotor with 
a finite number of blades. For example, with the optimum distributions 
of circulation over infinitely-bladed rotor disks derived from the 
general momentum theory, the effect of a finite number of blades is 
obtained approximately by using Prandtl ' s formula (Ref. 11). Also, 
Greenberg and Kaskel's flowfield results for the infinitely-bladed rotor 
has been used by Erickson (Ref. 6) in a prescribe-? wake analysis. The 
research reported in the present work utilizes the simplicity offered 
by the actuator disk concept, but incorporates certain effects of slip- 
stream rotation which have been neglected in previous theories for rotor 
performance . 
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III. GENERAL FORMULATION 


Consider a rotor hovering OGE in an incompressible and inviscid 
fluid. Following Prandtl's lifting line theory, each rotor blade is 
represented by a blade -bound vortex filament directed radially. If 
the strength of this vortex filament, i.e., the blade circulation, 
varies with the radial distance from the rotor axis, then there exists 
a vortex sheet trailing each blade. Within the context of the actuator 
disk concept, the total circulation of the several blades of the rotor 
is considered to be distributed uniformly over an infinite number of 
thin blades, each with a vanishingly small chord width. The discrete 
vortex sheets trailing the lifting lines (blades) are replaced by a 
volume distribution of vorticities within the slipstream. The problem 
is then time -independent and axi-symmetric about the rotor axis in a 
reference frame at rest relative to the fluid at infinity. That is, 
in a cylindrical coordinate system (r, 0, z) with the origin at the 
center of the disk and the z-axis coinciding with the rotor axis end 
pointing downstream, the flew parameters are all independent of the 
9-coordinate and of the time t. The tangential, or 0-, component of 
the velocity vector, however, is not zero within the slipstream. 

The general flow features for the hovering rotor OGE are shown 
in Fig. 1. For the present problem, the entire flow is induced by the 
rotor. Consequently the entire region outside of the slipstream and 
the blade row is "upstream" of the rotor. 
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The differential equations describing the flow are the familiar 


continuity and momentum equations : 


( 1 ) 


and 


(q • 7) q = - - 7 p + F 
P 


( 2 ) 


where F stands for the equivalent body force exerted on the fluid by 
the blade elements. The boundary conditions are '.hat q vanishes infinite- 
ly upstream of the rotor and is independent of z infinitely downstream, 
i.e., in the ultimate wake. 

For axi-symmetric flows, Eq. (l) can be rewritten as 


l ± (ru ) + & * 0 
r fcr v dz 


( 3 ) 


where u and w are the radial and axial velocity components respectively. 
Therefore a Stokes' stream function t|t(r, z) exists such that 


and 


1 & 

u rr~ 

r 3z 


„ I$± 

r 8r 


( 4 ) 


( 5 ) 


The flowfield is therefore completely determinate if f and the tangential 
velocity component v are known. 

Introducing the vorticity vector <2 defined by 


(a » 7 x q 
2k 


( 6 ) 



( 7 ) 




one obtains the components of oj in terms of ijr and v: 


* dv 
?a - 5 ? 




( 8 ) 


and 


„ 1 ^ t \ 

C*?& (rV> 


(9) 


The momentum equation can be rewritten in terms of u) and the 
total head H as 


q x w = ?H - F 


( 10 ) 


where 

H - £ + \ (11) 

P d 

Outside the blade row, F = 0 and therefore q • vH = 0 and uJ • vH = 0. 
Consequently, the gradient of H is perpendicular to both the velocity 
vector and the vorticity vector. In other words, the total head of the 
fluid is a constant along each stream tube except at the blade row 
where it may undergo an abrupt change. Upstream of the rotor, there- 
fore, the total head is everywhere equal to the static pressure of the 
fluid far upstream divided by the density. Inside the slipstream, the 
total head is a function of the stream function only; 

p 

H ■ -* upstream of the rotor (12) 
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f(ijr) inside the slipstream 


Noting further that outside the blade row, since F » 0 f*nd the 8- 
component of VH is zero because of axi-symmetry, one has 

(qx 5) • e e = o 

or 

q • 7 (rv) a 0 ( 13 ) 

Consequently, outside the blade row, the gradient of rv iB perpendicular 
to the velocity vector. In other words, the angular momentum of the 
fluid remains a constant along each stream tube except at the blade 
row where it may undergo an abrupt change. Since the angular momentum 
of the fluid is zero far upstream of the rotor, the angular momentum 
is zero everywhere outside the slipstream. Within the slipstream, the 
angular momentum is a function of the stream function only: 

vt ■ 0 upstream of the rotor (14) 

« g {if) inside the slipstream 

Consider the fluid passing the disk at r = r^ as shown in Fig. 1. 
Equation (l4) states that the tangential velocity component v of the 
fluid is zero on the upstream side of the disk and may be non-zero on 
the downstream side. Let v^ be the tangential velocity component of 
the fluid on the downstream side of the disk. There exists then a 
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discontinuity in v at the disk described by a distribution of radially 

directed vortices at the disk of strength v^. The circulation of the 

blade-bound vortices at r. is therefore 2nr,v.. The tangential veljcity 

d ad 

at the disk is in the principal value sense. 

Within the context of the actuator disk concept, the blade-bound 
vortices are radially directed and axi -symmetrical. Therefore, 
the axial and radial velocity components of the fluid are both continuous 
at the disk. The continuity of the axial velocity at the disk is of 
course consistent with the law of conservation of macs. 

Kutta- Joukowski theorem gives the lift (or thrust) and drag 
forces on the annular element of the disk ds^ = 2rrr d dr d as follows: 

dT = p (nr d - v d /2) v d ds d (15) 

dD = p w d v d ds d O- 6 ) 

The rate of work done, or power expended, by the disk element ds fl is 
therefore 


dP = 0r d dD 

= ^ w d V d r d ^d (17) 

The total thrust and total power of the rotor are obtainable by inte- 
grating Eq.8, (15) and (17) respectively over the disk. 

The mass flow rate through ds d is piw d ds d . Therefore, from Eq. 
(17), one finds the increase in total head across the rotor disk to be 

ftv. r.. The total head on the downstream side of the disk is therefore 
d d 
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+ n v d r d 


(18) 


H, 


pap 

P 


According to Eqs. (12) and (14), H and vr are both functions of f only- 
in the slipstream. One therefore has, from Eq. (l8), in the slipstream 


dH « d / N 

a* ■ n i* (vr) 


(19) 


The radial component of Eq. (10 ) gives 


„ _ 3H dH a* 

V C - wTl - ^ £ 

By using Eqs. (7), (9)> (19 ) j and (20) - one obtains 


( 20 ) 


11 s ■ (nr - v) 


( 21 ) 


Consider now Eq. (8) which relates T] to i|r. This equation is 
readily identified with the problem of circular vortices having the 
z-axis as a common axis. In many classical treatises on hydrodynamics 
(e.g., Ref. 25), the special case of an axi-symmetric rotational flow 
with v, 5 , and £ all zero is discussed. The present problem is more 
general in that v, % and £ are non-zero within the slipstream and, 
furthermore, v(=> £ v^ ) and 5 (blade-bound vortices) are non-zero at 
the rotor disk. The differential equations relating u, w, and T) to % 
for the special case of zero v, ?, and £, however, are identical to Eqs. 
(4), (5), and (8) derived for the general problem. This fact is not 
surprising since for the general axi-symmetric flow the continuity 
equation does not contain v and neither the radial vorticity distribution 
nor the axial vorticity distribution induces a velocity in a meridional 
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plane. As a consequence of this fact, certain kinematic relations 
between u. w, T] and i|r available in the classical treatises for the 
special case of zero v, 5, and Q are directly applicable to the present, 
more general, problem. In particular, the differential equation (8) for 
i|r may be immediately recast into an integral representation for \jt : 

f(r, z) = JJ G(r, z; r', z') T]' dr* dz* (22) 

D 

where the region of integration D covers the entire slipstream, G is 
the stream function at a field point (r, z) associated wiih a circular 
vortex filament with unit strerr ;h of radius 1 * and centers about the 
z axis in the plane z « z'. T)' is the vorticity vai.ue in the slipstream 

at r ' , z ' . 

Several alternative forma of this Green's function are available 
in the literature (e.g.. Ref. 25) and different authors developed 
different forms in their studies of the rotor and propeller problems. 

For example, an integral form involving Bessel functions of the first 
kind was utilized in Ref. 3 and a Legendre function of the second kind 
and the degree \ was utilized in Refs. 4 and 5. A convenient form to 
use is one involving complete elliptical integrals (Ref. 25, 26): 

_ — - - — - - 2 

a \ (r + T'f + (z - z') 2 { ( 1 - \ ) K(k) - E(k) } (23) 


where 


[ 


4rr' 


(r + r* y + (z - z' ) 


4 

2 j 


(24) 



1 


(25) 


K(k) = J 


h/2 


° V 1 - k 2 sin 2 a 


da 


»tt/2 j p p 

E(k) = J V 1 - k sin a da 


( 26 ) 


The functions K and E axe complete elliptic integrals of the first and 
the second kind respectively. The properties and numerical values of 
these functions are readily available in mathematical handbooks . 

The integral, representation for i|r, Eq. (22) is a kinematic 
relation allowing t)r to be determined from any known distribution of 
T\ in the slipstream. It replaces the differential Eq. (8) together 
with certain boundary conditions. By putting Eq. (21 ) into Eq. (22), 
one obtfins an integral representation equivalent to the ones derived 
by T. Y. Wu and by Greenberg, et. al. (Refs. 3> and 5) using Green's 
functions. Greenberg et al. recognized the Green's function as the 
stream function due to a ring vortex. They subsequently used a vortex 
representation of the slipstream and developed an efficient method for 
computing the flowfields induced by the rotor. 

During the course of the present investigation, several closed- 
form analytical, solutions have been developed for the stream function 
and for velocities induced by several uniform and linear distributions 
of circular vortices over finite circular cylinders and annular surfaces. 
The solutions are expressed in terms of complete elliptic integrals 
and are presented in Appendix A of this report. 


) 
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The availability of these analytical expressions contributes to the 
efficient and accurate computation of axi- symmetric flows. 
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IV. OPTIMUM PERFORMANCE THEORY 


The Ultimate Wake 

Consider the rotor hovering OGE. As shown in Fig. 1, the slip- 
stream contracts downstream of the rotor and forms the ultimate wake 
far downstream where it no longer contracts. Thus, in the ultimate wake 
the radial velocity u vanishes and the flow is independent of the axial 
coordinate z. The fluid passing through the disk element ds^ passes 
through the annular element ds^ = 2rrr w dr w in the ultimate wake as 
shown in Fig. 1. By the law of mass conservation, one has 


Equation (l4) gives 


w. ds. = w ds 
d d w w 


( 27 ) 


v. r. = v r 
d d w w 


The radial component of Eq. (2) gives in the ultimate wake 



( 28 ) 


( 29 ) 


This equation states simply that the centrifugal force of the fluid is 
balanced by the radial pressure gradient in the ultimate wake. Consider 
a spherical control volume of radius R centered at the origin of the 
coordinate system. The control surface consists of two parts and Sg, 
where S 1 is a cross section of the slipstream and Sg is the remaining 
portion of the sphere. At large distances from the rotor disk, the 
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action of the rotor on the fluid outside the slipstream is equivalent 
to that of a sink, the strength of which is equal to the volume rate 
o; flow passing through the disk, and is finite. Therefore, in the 
limit as ft -• ®, p = p^ + 0(l/ft^). Therefore the static pressure at the 
rim of the ultimate wake is p . Equation (29) therefore gives 

JR v 2 

p. - = m w Hr (30) 

r w 


Equation (12) requires H w = Hj. It follows from Eqs. (l4), (18), and 


(30) that 




P _R v 2 

= — + fl v r +r w — dr 
p w w J r 


w 


(31) 


However, in the ultimate wake, since u^ = 0, one has 


2 2 

T) V W 

H _ + Jl_ + JL. 

\ “ D + ~ + ~ 


(32) 


Therefore 


2 2 D 2 
W V J* V 

w _ w . r w w , 

-7T— = n v r — =— + dr 

2 w w 2 J r 


w 


At the rim of the ult mate wake, 


(33) 


w_. 


w o n w o 

“ = n v w R w-"T 


(34) 


Consider again the control surface described earlier. The momentum 
theorem gives the following expression for the total thrust of the rotor 
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T = J (p + pw 2 ) n • e z dS 1 + J (p + pt 2 ) n • dig {&) 
S 1 S 2 

where n is the unit outward normal vector on and Sg. As R -> ®, t: 
integrand in the second integral goes to p^ + 0(1/1^). Therefore the 
second integral give** - _ co (ttR w ). Equation (35) therefore gives, with 
the use of Eq. (30 ) and Eq. (33), 

T - ip L [ ' r w (2(lr w - V + \ 2 ] (36) 

S 1 

The total power expended is, by using Eqs. (17), (27), and (28), 

P - & L w w T w r » (37) 

S 1 

With a known radial distribution of the tangential velocity v 
in the ultimate wake, Eq. (33) establishes the distribution of w^. 
Equations (36) and (37) then give the thrust and power of the rotor. 

A perturbation in the tangential velocity distribution therefore alters 
the thrust and the power. In the following section, a criterion for 
optimum performance is established by analyzing the changes in power 
and thru t resulting from the perturbation of the tangential velocity 
distribution. This criterion is first given in terms of the velocity 
distributions in the ultimate wake. These velocity distributions will 
subsequently be related to the optimum distribution of bound vorticity, 
thrust, and power. 
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Criterion for Optimum Performance 

An optimum rotor is a rotor giving a specified amount of thrust 
while expending a minimum amount of power. Consider a given distribution 
of tangential velocity v w in the ultimate wake. If a perturbation Av^. 
in the tangential velocity is introduced as follows: 


Av = 
w 


0 , for 0 < r w < a 

, for a < r w < a + Ar w 


0 , for a + Ar < r < R 

’ w w w 


(38) 


then, according to E\. (33 )> the resulting change in the axial velocity 
distribution is 


Aw = < 
w ^ 


2v A r e/w , for 0 < r < a 

wa w ' w w 

(na - v^Jae/w^, for a < r^ < a + Ar. 


w 


w 


w 


, for a + AF W < r w < 


(39) 


where the subscript a indicates the value of v is evaluated at r = a. 

w w 

Only terms of the lowest order in ea/v and in Ar / a are kept in each 

wa n 

interval of interest. 

Equations (36), (37) > together with (38) and (39) gives the 
changes in thrust and power: 


and 


AT = 2rTp(eAr w ) a (2fla - v^) 


AP » 2TTffi(eAr w ) {a- 3 v wa (r>a - v^J/w^ 


(40) 
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1 


(41) 


a^w 


wa 


+ 2 v 


wa 


/ 


,a v 


w w_ 

w 

w 


dr } 

•r.r J 


w 


Equations (4o) and (4l) give 


AP = h(a) AT 


(42) 


where 


h(a) = 


0 fav (fla - v )/w + aw + 

L wa wa ' wa wa 



2fla - v. 


wa 


(43) 


is a function of a, the location where the perturbation of tangential 
velocity is introduced. 

Suppose that a perturbation in tangential velocity is introduced 
at a 1 with resulting changes in power and thrust, t£' 1 and A^. Suppose 
a second perturbation in tangential velocity is introduced at a^ with 
resulting changes in power and thrust, AP 2 and AT 2 * suitably select- 
ing the values of eAr^, at a^ and a^ , the total change in thrust due to 
the two perturbations is made to vanish. Then ATg = - AT^ and therefore 

AP = AP X + AP 2 = [h(a 1 ) - h(a 2 )] (44) 


Since AT-^ can be made either positive or negative by choosing the sign 
of e in Eq. (4o), it will be possible to obtain a negative value of AP 
if h(a^) ^ h (a 2 ). In other words, it is possible, by altering the 
tangential velocity distribution, to develop a given amount of thrust 
while expending a smaller amount of power. The tangential velocity 
distribution v^ under consideration therefore does not give optimum 
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performance . 

It follows therefore that the optimum distribution of tangential 
velocity v w is one which gives h(a^) = h(ag). Since and ag are 
arbitrarily selected radial locations, the criterion for optimum per- 
formance is 


h(r w ) = constant 


(45) 


In other words, 
r 


r Wv p 

_W 2 _ 

r d 
w 

‘ JL ( 

■**» . \ 

w 

1 - _W . 

Or - v 
w w | 

J — r ar - 

o w 
w 

2 

_ r ' 

w 


' V 

w 

w J 

w 


(46) 


where N is a constant. 

Computation Procedure 

Equation (33) relates w w to v . As a consequence, Eqs. (33) and 

(46) permit the evaluation of the optimum distribution of the tangential 

velocity v in the ultimate wake for specified values of N and R . 

w « 

Alternatively, Eqs. (33) and (46) are re-expressed in non-dimensional 
forms : 


~ 2 ~ 2 . ~ 2 

w v J. v 

v ~ — w r w 

~ = V r - dr 

2 w w 2 ~ 


(47) 


w 


and 


w v 

J. f 


~ 

r dr 




2 r 


r - v 


LTIf -i. i \ _ _w _ v 

/ ~ ~ 


w 


w 


w 

w 


- ] (M) 


where 


w 

n*w 
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i: 


w 


w 


W fJR 


w 


v R 


w 


(49) 


and 


». J- 


From Eq. (4?), one obtains at * 1 


* 


w = [ V (2 - V )] 2 
w L w w /J 

o o o 


(50) 


where the subscript "o" denotes the velocity components at the rim of 
the ultimate wake. Differentiating Eq. (47) with respect to r w yields 


dw dv 

— £„ (? - y )(v + r — * )/r w 

or dr 

w w 


(51) 


From Eq. (48) and (50), one obtains at r 


w 


(3 ■ 2 v v )[v w / (2 - v w )]^ + 2v w J 1 (v w r 2 /^) dr 


Sf 


o o 


o o 


2 - v 


(52) 


w 


Differentiating Eq. (48) with respect to r^ and using Eq. (51) yields 


dv v [3w ^ - 2fJ(3r v )w ^/r + v (5r - 2v )w - v 2 (r - v ) 2 ] 

w w w w w w' V ' W W v W W / W W v W w' J ,_„v 

“5 (53) 

j ^ r— rrr /m*t -j rr f r T i :t it * 

^w 


_W W L - w ~ w w 

r [w k -2t KT 3 - v (r - 2v ) w * + v 2 (r - v )*] 
W L W W W W V w w w w J 
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If the value of is specified, then w is determinate from 
o o 

Eq. (50). If, further, the value of fif is given, Equations (51) and 
(53) can be re-written as a system of two first order ordinary differ- 
ential equations with "initial" values specified at the point = 1, 
i.e., the rim of the ultimate wake. Various numerical methods are 
available for the solution of these equations. In the present work, 
a fourth order Runge-Kutta method is used to determine the distributions 
v and w frcm Equations (51) and (53) for specified values of v and 

w w w 

~ 0 

!Rf. An iterative procedure is used with which a value of N is estimated 
for the initial iteration. The distributions v^ and w^., obtained using 
this value of fT from Eqs. (51) and (53) are placed into Eq. (52) to 
compute a new value of ft for the second iteration. The iterations con- 
tinue until the value of fl changes insignificantly between two successive 
iterations. In this manner, the optimum ultimate wake flow corresponding 

to a specified value of v is established. 

w 

o 

Contraction Ratio and Performance 

The figure of merit of a rotor is defined by 


1 V 


J? C P 


where C T and Cp are respectively the thrust coefficient and the power 


coefficient defined by 


C T = T/TTpfi 2 R d 4 


C p = P/npfi 3 R d 5 


29 




For a given radius of the ultimate wake, R^, and a given angular 

velocity fl of the rotor, the computed optimum distributions of v and w 

w w 

together with Eqs. (36) and (37) yield the corresponding thrust and 
power of the rotor. In order to establish the performance of the rotor, 
i.e., C^, Cp and M, however, the disk radius must be determined. 

In theories of propeller performance, it is customary to write 
the thrust element in terms of the wake velocities as follows: 


dT = £ P[(2flr w - v w )v w + w y d ] ds w 


The acceptance of Eq. (57) permits the radius of a given stream tube 
on the disk, r^, to be related to the radius of that stream tube in 
the ultimate wake, r . Equating the right sides of Eqs. (15) and (57) 
and using Eq. (28) gives 


r (2nr - v )v + w ] r, 
1 v w w 7 w w J c 

2 

(2flr . - v r )v 

' d w w' w 


The boundary condition for this differential equation is 


r , » 0 at r = 0 
d w 


With known distributions v w and w^, the quadrature of Eq. (58) gives 
as a function of r^ and thus, for the given value of R w , the value of 


Equation (57) obviously satisfies Eq. (36). However, it is not 


the only expression that satisfies Eq. (36). Any function b(r ) which 

Jt w 

w 

gives j bdr w = 0 can be added to the right side of Eq. (57) and the 
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resulting expression for df still satisfies Eq. (36). It has been 
shown in Ref. 29, in fact, that no solution of Eq. ( 58 ) exists which 

satisfies the boundary condition (59). Accordingly, in the present 
research, an approach which does not utilize Eq. ( 58 ) is developed. 
In this approach, the contraction ration 


0 ■ V r d 


(60) 


is taken to be a constant independent of the radial location (r w or r^). 

In general, the contraction ratio is a function of the radial 
location. The results of Ref. 2, however, indicate that no serious 
error is introduced by taking the contraction ratio to be a constant 
for heavily loaded free-running propellers. Thus it is expected that 
the assumption of constant contraction ratio is reasonable for the 
hovering rotor. 

By taking c in Eq. ( 60 ) as a constant, Eqs. (15) and (28) yield 


Integrating Eq. 
equation to the 
terms, 


/ Or v- 

dT - p ( -r - v ) v . < 61 > 

c 

( 61 ) and equating the right siae of the resulting 
right side of Eq. ( 36 ) gives, after rearrangement of 


c 


[ 


1 + 




r ds 
w w 


( 62 ) 


Noting that, using Integration by parts, 
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% 

r 2 

w r dr 
J v w w 


2_ 2 

w R 
w o w 


R . 2 

w dw_ 


l f" w 2 

' 4 J. £T r w dr « 


O V 


( 63 ) 


Differentiating Eq. (33) gives 


2 2 2 
dw.. dv.. dv v 

t 3T “ + ^ “ - i ^ — “ 


w *^w dr 


dr 


w 


(64) 


Placing Eq. (64 ) into (63 ) and noting that 


w dv 


f r-^r 3 dr = v R 3 - 3 [ v r 2 dr 
J dr., w w ww J ww w 


(65) 


and 


% 2 

J — r w dr . 

o w 


2 2 

v R - 2 
w w 


lAr 2 
J v w r w 


dr 


( 66 ) 


One then has 


p w 2 


w r dr = 


w 2 R 
w w 

0 


J n a v** — A 

.www 2 


- <v„ 3 * ± + 20 r vvH < 67 ) 


Using Eq. (34) then gives 


P w 2 I* w 2 

J w » r w dr w* 20 J v w r w 


( 68 ) 


Consequently, Eq. (62) becomes 


c ■ 


R 


w 



(69) 
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It is noted that this contraction ratio is precisely the value predicted 
by the axial momentum theorem. 

With c known, R^ immediately follows from the specified value 
of R w (R^ « J~2 R w ). The axial inflow velocity and the circulation at 
the rotor disk are related to the axial and tangential velocities in 
the ultimate wake. The thrust and power coefficients as well as the 
figure of merit can be easily evaluated from the known flow conditions 
in the ultimate wake. In particular, once an optimum distribution of 
v is established, the corresponding distribution of w^, w^, and P, the 
values of C_, C , and M immediately follow. 

I X 

The optimum rotor performance theory described here has been 
generalized to the case of a rotor (or a propeller) in axial motion. 

The analysis for this general case is presented in Appendix B of this 
report . 
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/. R0 T R INDUCED FLOWFIELD 


A number of methods have been developed by various authors for 
the solution of Eqs. (21 ) and (22), with specified distribution of 
disk-bound circulation (Refs. 3> 4, 5» 23, 24). The approach ch<~ an 
for the present work is the numerical method described by Greenberg 
et al. (Refs. 4 and 5). 

For the hovering rotor problem, the domain of integration in 
Eq. (22) is the entire slipstream and is not known a priori. In fact, 
the shape of the stream tubes bounding and within the slipstream is 
to be established as a part of the solution. Following the approach 
of Greenberg and Powers (Ref. 5)j the specified distribution of circula- 
tion at the disk is approximated by a piecewise uniform distribution 
as shown in solid lines in Fig. 2. The disk radius is divided into L 
segments at radii r^ = R q , R^, . . . R^, with R q = 0, and R^ = R^. In 
the l th segment, R < r , < R., the circulation is taken to be a 

-L Cl 4* 

constant value F^. The trailing vortices are considered to be shed 
only at the radial locations R^. The slipstream vortex system is then 
represented by L discrete vortex tubes . The <tth vortex trails the disk 
at r^ r r^. Equation (22) becomes 


t( r > 




(70) 


where t_£ is the radius of the tth vortex tube and y^' is the strength 
of the -tbh vortex tube, t ' and y ' are both functions of z'. 

v 4 * 
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With the above representation of the slipstream vortex system, 
the vorticity in the regions between the vortex tubes is zero. The 
angular momentum vr of the fluid is constant between the vortex tubes 
and changes discontinuously across the tubes by the amount (r^ +1 - 
/2tt. Consider for the moment a continuous, though rapid, change of 
the angular momentum in the interval r = t^tor = i^ + £r. The 
corresponding volume distribution of the tangential vorticity TL is 

V 

given by Eq. (21 ), i.e. 

\ = - «* - v) (71) 

As At approaches zero, one obtains a vortex tube with strength 

r 


The angular momentum vr at r = t^ is, in the principal value sense, 
equal to \ [(vr)^ + (vr)^ j , or equivalently (r^ + F^ +1 )/4 tt. Thus 
Eq. (72) may be re-written as 



Using Eq. (5), one then obtains 


r - r 

y l l 2 tt 


[»- 


4rrt 




+ r 


i+i 


)] 


(74) 


where w. is the axial velocity on the tth vortex tube in the principal 
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value sense. 

Equation (74) has been derived in Ref. 5 using the Bernoulli’s 
equation and the condition that the vorticity in the slipstream convects 
with the fluid motion. The present, alternative, derivation of this 
equation points out that, instead of considering the trailing vortex 
tubes to originate at the disk radii R^, it would also be appropriate 
to lump the vortices shed within each interval R, < r . < R, into a 
vortex tube originating at the mid-point of that interval, i.e., at 


r d = ^(R^^ + R ^)* This mid-point representation is shown in dashed 
lines in Fig. 2. While the end-point representation used in Ref. 5 
describes the disk bound circulation well, it appears that the mid -point 
representation would describe the slipstream vortex system more accurately. 
With the mid-point representation, if the circulation goes to zero contin- 
uously at r^ = R^, then the slipstream boundary is no longer identified 
by the outermost vortex tube. For this situation, the mid-point repre- 
sentation must be modified for the outermost interval if the slipstream 
boundary is to be identified. The mid -point representation also presents 
seme difficulties for the innermost interval. In Ref. 5, attempts were 
made to provide a cut-out, i.e., to let = 0, representing the finite 
hub of the rotor. These attempts failed to produce convergent results. 

The use of the mid-point representation for the innermost interval is 
computationally equivalent to the use of an end-point representation and 
assign = 0 for the interval 0 sr^ s R^/2. Thus it is expected that 
the use of the mid-point, xepresentation for the innermost interval would 
also present convergence problems. 


/ 
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A possible method of removing the above noted difficulties is to 
use continuous distributions of blade-bound vortices, or circulation, 
for the innermost and the outermost intervals. It should be noted 
that in Refs. 4 and 5, it was pointed out that with the end point repre- 
sentation, the strength of the outermost vortex tube becomes singular 
at the tip of the slipstream. Furthermore, if the circulation for 
the innermost interval is non-zero, then the tangential velocity v 
becomes infinitely large at r = 0 and the rotor's thrust diverges. 
Therefore there existed also difficulties in using end-point representa- 
tions for the innermost and outermost intervals. Again, a possible 
method of removing these difficulties is to use continuous distributions 
of blade-bound vortices for the innermost and outermost intervals. 
Inasmuch as satisfactory results were obtained by the use of the method 
described in Ref. 5 for the major portion of the flowfield of concern, 
no attempt has been made to incorporate the "continuous" representation 
of blade-bound vortices. 

Since the vortex tubes in the slipstream coincide with the 
stream tubes, the application of Eq. (70) at each vortex tube m gives 

♦m=, s I 0(t m , t y’ 4*’ (75) 

«Cr= 1 O 


where is a constant. Equations (5), (70 ) and (74) give 


L 

E 

- 6-1 



.. + » 


z*) 


dz’ = 


t O' 
m v m 




2tt V 
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[«- 


4nt_ 


< r n + 


m 



where 


G r (r, z; r', z' ) 1 { K(k) 

2rrr ' J (z’-z) 2 + (r’+r) 2 



1 + 


2r (r-r* ) 
(z’-z) 2 + (r’+r) 


o ] E(k) } 


(77) 


Equations (75) and (76) constitvce a set of 2L coupled integral 
equations containing the 2L unknown functions t^(z) and y^( z ) . An 
iterative method of solution of this sec of integral equations is 
described in Ref. 5. 
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VI. RESULTS AND DISCUSSION 


Optimum Performance 
a. Optimum Circulation Distribution 

Numerical results for the optimum radial distribution of circu- 
lation and inflow velocities at the rotor disk and the associated per- 
formance parameters have been obtained for l6 values of the thrust 
coefficient between 0.001 and 0.050. This range of thrus* coefficient 
was considered to bracket the normal operating range of hovering rotors. 
In Ref. 27, the optimum radial distributir is of the circulation and of 
the axial inflow velocity at the disk axe presented for the 16 values 
of the thrust coefficient. It is noted in Ref. 27 that for small values 
of the thrust coefficient, the optimum circulation and tne axial inflow 
velocity are nearly independent of the radial position except neax the 
axis. For higher values of the thrust coefficient, the deviations of 
the optimum circulation and of the axial inflow velocity from constant 
values extend over laxger regions and the magnitude of the deviations 
axe larger. 

The optimum distribution of circulation is presented for the 
thrust coefficients of 0.001, 0.005, 0.010, and 0.050 in Fig. 3- In 
Ref. 10, an approximate solution for the hovering rotor problem is 
given based on the general momentum theory neglecting the radial pres- 
sure gradient in the umtimate wake. The approximate distributions of 
circulation of Ref. 10 are computed and compared with the present, more 
exact, solution in Figure 3. The comparison shows that 



the approximate solution of Ref. 10 deviates significantly "Prom the 
present result, particularly in the inboard region of the disk. The 
deviations are more severe for larger values of the thrust coefficient. 

b. Axial Inflow Velocity 

In Fig. 4, the optimum distributions of axial inflow velocity 
obtained in the present study are presented and compared with the approxi- 
mate solution of Ref. 10 for several values of C T . It is noted that 
the approximate solution of Ref. 10 gives a zero value of the axial 
inflow velocity at the axis of the disk. It has been shown in Refs. 

27 and 29 that this zero axial inflow velocity is incorrect. 

c. Figure of Merit 

Insofar as the overall performance of the rotor is concerned, 
the approximate method of Ref. 10 yields a figure of merit more 
optimistic than the more exact results of the present theory. Figure 
5 shows that the computed figure of merit based on the present theory 
is lower than that baaed on the approximate theory of Ref. 10, which 
is in turn much lower than the ideal figure of merit of one. The 
deviations between the results increases with the value of the thrust 
coefficient . 

d. Effect of Finite Number of Blades 

While the present investigation is primarily concerned with an 
infinitely bladed rotor, the effect of a finite number of blades may be 
determined, for example, by extending Lerbs' work (Ref. 17) on heavily 
loaded, free running propellers. 

An estimate of the effect of number of blades on the optimum 
distribution of circulation is provided by modifying Prandtl's approxi- 
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mate method (Ref. 11 ) of calculating the "tip loss" factor and applying 
the result to the hovering rotor case. The procedure is outlined in 
Ref. 27. The results are shown in Figure 6 for rotors with 2, 4, 
and 6 blades at = 0.010. Based on these results, the corrected 
thrust coefficients for 2, 4, and 6 blades are respectively 0.00759, 
0.00865, and .00906. The corresponding corrected power coefficients 
are .000547, 0.000624, and 0.000655. The figures of merit are 0.855, 
0.911, and O.931. In comparison, without the tip loss correction, the 
power coefficients for the infinitely bladed rotor are respectively 
(at C T = 0.00759, 0.00865, and 0.00906) 0.000479, 0.000585, and 0.000627. 
The figures of merits are 0.975, 0.973, and 0.972. 

Rotor Induced Flow 

Following the analyses described in the previous sections, a 
computer program was prepared to compute the flowfields associated 
with the optimum distribution of circulation. In these computations, 
emphasis was placed on the determination of the shapes of the stream 
tubes in the slipstream. As discussed earlier, the vortex tubes coin- 
cide with the stream tubes in the slipstream. 

The continuous circulation distribution was represented by a 
number of piecewise uniform circulation intervals and the trailing 
vortex tubes are considered to originate from the end-points of each 
interval as was done in Ref. 5. The value of the circulation within 
the iVfch interval was taken to be 


r 


1 




V T a V'd 


(78) 
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It shall be shown later that the thrust developed in the tth interval 
due to the constant circulation given by Eq. ( 78 ) is approximately 
equal to that due to the continuous distribution of circulation r. = 

2 " "a r d‘ 

a . Effects of Number of Intervals Used 

Three sets of computations were made for the case of an optimum 
rotor hovering OGE with a thrust coefficient of 0.010 using h, 5 } and 
7 constant circulation intervals. For the set of computations using 
4 intervals, the values of R^'s were assigned to be 0.00 R, 0.15 R, 

0.30 R, 0.60 R, and 1.00 R. For the 5- interval computations, additional 
R^ values were assigned at 0.80 R. For the 7-interval computations, 
additional R^ values were assigned at 0.4-5 R, 0.75 R> and 0.90 R. A 
comparison of the computed stream tube shapes based on 4-, 5-, and 7- 
interval representations is shown in Fig. 7. The comparison revealed 
that the three sets of computations gave practically identical stream 
tube shapes . The maximum deviations in the computed radius of the 
slipstream boundary was less than 0.2$ between the three sets of results. 
A similar comparison of the axial inflow velocities revealed that, 
except in the immediate vicinity of the additional R^'s, the results 
are virtually independent of the number of intervals used. These com- 
parisons suggest that a 4- or 5- interval representation is adequate 
for routine computations. 

b. Effects of Second Order Terms 

According to Eq. (15), if is such that the thrust developed 
in the tth interval due to is equal to tliat due to the continuous 

distribution of circulation f = 2tt v . r., then 
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The second terra on the left side of Eq. (79) represents a second 

order effect of the slipstream rotation. If this term is neglected, 

then one obtains the expression for r. given by Eq. (78). In general, 

2 

the maximum value of is approximately equal to 2nOR^ C^,. For usual 

applications, since C^, is much smaller than unity, the second term on 

the left side of Eq. (79) is negligible for all intervals 4/1. For 

the first interval, 4=1, one has R , , = R =0 and the second term 

becomes infinitely large. Thus, Eq. (79) is not appropriate for the 

first interval. It should be noted that in reality, the area represented 

2 

by the first interval is equal to ttR^ and is a small fraction of the 
total disk area ttR^ . For example, by taking R^ = 0.15 R^, the area 
represented by the first interval is 2.25$ of the total disk area. 

Since the value of the circulation in the first interval is also small, 
the thrust developed in the first interval is only a s ma ll fraction of 
the total rotor thrust. The piecewise constant circulation as deter- 
mined by Eq. ( 78 ) therefore gives a total rotor thrust very nearly 
equal to the total rotor thrust due to the continuous circulation dis- 
tribution T = 2tt v . r. . 

d d 

The second term inside the brackets in Eq. (74) also represents 
a second order effect of the slipstream rotation. This term again is 
negligible for all except the innermost interval. In this investigation , 
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the role cf this second order term in defining the vortex tube shapes 
is examined by computing, for the C T = 0.01 case using 7-intervals , 
the vortex tube shapes neglecting and accounting for this second-order 
effect. The results, shown in Figure 8, indicate that the neglect of 
this second-order effect leads to a slightly smaller ult ima te wake. 

The effect is small for all except the innermost streamtube. 
c. Restriction on the Radius of the Innermost Tube 

Equation (79) shows that if as r^ -* 0, the continuous distribution 
of circulation, r, goes to zero as r“, with a> 0, then the thrust 
developed in the first interval is finite. If the continuous distri- 
bution of circulation is replaced by a constant, non- zero, circulation 
in the first interval, then as r^ -• 0 the tangential velocity component 
goes to infinity as l/r^ and the thrust becomes infinitely large. 
Greenberg et. al. made several attempts (Ref. 5) to remove this diffi- 
culty by letting 1^ = 0. They pointed out that this is equivalent to 
providing a cut-out to represent a finite hub of the rotor. These 
attempts, however, failed to produce any convergence solution. In the 
present investigation, it was observed that the convergence of the 
numerical solution was also sensitive to the value of R^ used in the 
computation procedure. Convergence was not obtained with R^ = 0.10. 
However, with R^ => 0.15, convergence was obtained. The difficulties 
experienced is attributable to the fact that the innermost vortex tube 
was situated in an essentially ciead-air or recirculating flow region, 
as can be shown by examining the axial velocity in the ultimate wake. 

Since the vorticity between two adjacent vortex tubes is zero 
and since the radial velocity is zero in the ultimate wake, the axial 



velocity in the ultimate wake is constant between two vortex tubes. 
Let W. be the constant ultimate axial velocity in the -tth interval 
(i.e., between the (-t-l)th and the -tth vortex tubes), then, in the 
ultimate wake 




W. - W 


t+1 


(80) 


^ ^ + 


(81) 


Applying Eq. (7k) in the ultimate wake then yields 

= V l+1 + 2 ‘ r *l> " 1+TT 2 t 2 ‘ T l+1 ) ^ 

Noting that the axial velocity and the circulation is zero out- 
side the ultimate wake and neglecting the second order terms in Eq. (82) 
gives 

- S r t 

and in particular 


® a Q r 

1 TT r i 


Consequently, if one lets ■ 0 as was attempted in Ref. 5, 
then a o and the net flow passing within the innermost stream tube 
is zero. Thus the innermost vortex tube contains a dead-air or recir- 
culating flow region. Such regions obviously are not amenable to the 
iterative proceedure used here. 


It can be shewn that the inclusion of the second order terms in 



Eq. (82) increases the values of W. Therefore, is non-zero even if 
one lets « 0. This non-zero value of W^, however, is small and the 
innermost stream tube still contains essentially a dead-air region. 
Similarly, since the continuous distribution of circulation approaches 
zero near the axis of the rotor, the value of is small if the value 
of is small. As a result, becomes very small and the innermost 

stream tube again contains essentially a dead-air region. These 
observations were substantiated by the numerical results obtained in 
this investigation, 
d. Stream Tube Shapes 

In addition to the case of = 0.010, computations were made 
for the previously discussed optimum distributions of circulations at 
C T = 0.001, 0.005, 0.050. For these additional computations, the 4- 
interval representation was vised. The comjuted slipstream shapes are 
presented in Figs. 9 for cases of 0.001, 0.010, and 0.050. The shape 
for the case of C T = 0.005 lies between those for C T * 0.001 and 0.010 
and is not shown. The results show that with increasing values of C T , 
the "dead-air" region near the hub becomes more prominent . The contraction 
ratio is nearly independent of the radial location. For the cases studied, 
this ratio deviated from the universal value of 0.742 by less than one 
half of a percent. This value is somewhat larger than the value 0.707 
predicted by t..e momentum theory. The computed axial velocity in the 
ultimate wake, also shown in Figs. 9, for C T * 0.010, is in good agree- 
ment with that predicted by Eq. (82). Furthermore, for the smaller 
values of C T , the computed axial velocity is in good agreement with the 
optimum ultimate wake velocity. For the case of ■ 0.050, however, 
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the computed axial velocity is significantly lower than the optimum wake 
velocity near the axis of the rotor. 

In addition to the optimum circulation distribution cases, com- 
putations were made using the circulation distribution given in Fig. 6 
for the two-bladed rotor with tip-loss correction. This special distri- 
bution of circulation shall b'i referred to as off-optimum circulation since 
the present computation procedure is specifically designed for infinitely 
bladed rotors. 

The slipstream geometry for the off -optimum circulation case is 
compared with the optimum circulation case of C T * 0.010 in Fig. 10. 

The comparison shows that the ultimate wake radius is smaller for the 
off -optimum case than for the optimum case. 

The slipstream boundary shapes for a rotor with constant circu- 
lation hovering in ground effect at heights and 2R^ above the ground 
plane are computed and shown in Fig. 11. Also shown in Fig. 11 are 
the experimental results of Ref. 18 obtained from smoke pictures for 
a two-bladed rotor hovering at a height of R d above che ground plane. 

It appears that, within the limitations of the infinitely bladed model 
for the rotor, the computed results provide * reasonable description of 
■the wake boundary. Computations were also made for a rotor hovering at 
a height 10 R d above the ground plane. It was found for this case, 
the flowfield at distances 3R d or more above the ground plane is nearly 
identical to that of a rotor hovering OGE. 
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VII. CONCLUSIONS 


Studies of the optimum performance and of the induced potential 
flowfields of rotors and propellers axe made on the basis of a "rotating 
actuator disk" concept. By modifying and adopting certain principles and 
techniques of propeller aerodynamics, a new criterion for the optimum 
performance of rotors and propellers in axial flight is obtained. The 
analysis leading to this criterion is more complete than the previous 
analyses in that the present analysis fully accounts for the effects of 
slipstream rotation and allows for the existence of a radial pressure 
gradient in the ultimate wake. For flowfield computations, the radial 
distribution of circulation at the rotor disk is represented by a piece- 
wise uniform distribution of circulation. The slipstream is represented 
by a set of axi-symmetric vortex tubes shed from the disk at the locations 
of circulation discontinuities. Existing numerical methods are utilized 
to compute the shape and strength of the vortex tubes that satisfy the 
required kinetic and kinematic conditions for the slipstream. Numerical 
results are obtained for hovering rote rs (static propellers ) in and out 
of ground effect for cases of optimum and off-optimum circulation dis- 
tributions at various thrust levels. 

Several conclusions of the present study are summarized below: 

1. Previous investigators comparing experimental data with 
theoretical results noted that the earlier theories, in which the effect 
of slipstream rotation is partly or totally neglected, generally over- 
estimate the optimum performance of hovering rotors in the higher thrust 
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coefficient range of practical interest. The present theory provides 
improvements in theoretical predictions by fully accounting for the 
effects of slipstream rotation. These improvements, which are not 
important in the case of a lightly loaded propeller, have significant 
influence on the optimum distribution of blade-bound circulation and 
hence power and thrust requirements of hovering rotors. The inaccuracies 
of the earlier theoretical results that are attributable to the neglected 
slipstream rotation effects increase with increasing thrust coefficient. 

Thus the improvements of the present theory are particularly important 
in the higher thrust coefficient range. 

2. The rotating actuator disk concept permits the calculation, 
in a straightforward manner, of the optimum distribution of circulation 
and inflow velocity over the rotor disk for the infinitely bladed case. 

The potential flowfield induced by the infinitely bladed rotor and the 
effect of finite blade number on the circulation distribution are 
estimated, by using available methods, for the optimum circulation dis- 
tributions . 

3. In the flowfield computation, the continuous circulation dis- 
tribution was represented by a number of piecewise uniform circulation 
intervals. The computed results show that representation of the continuous 
circulation distribution by as few as four piecewise uniform intervals 

is adequate for routine computations of the slipstream and of the stream- 
tubes within it. 

4. The computed slipstream shapes for rotors hovering OGE show 

that the dead-air region near the hub becomes more prominent with increasing 
value of the thrust coefficient. The slipstream contraction ratio is 
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nearly independent of the radial location and of the thrust coefficient. 
The computed value of the contraction ratio is somewhat larger than the 
value of 0.707 predicted by the axial momentum theory. 

5. Computations for rotors hovering in ground effect show that 
within the limitations of the infinitely bladed model, the computed results 

provide a reasonable description of the wake boundary. For a rotor 
havering at a height of ten rotor radii above the ground plane, the 
flowfield above about 3 radii distance from the ground plane is nearly 
identical to that of a rotor hovering 0(25. 

6. A number of explicit formulae useful in computing rotor and 
propeller induced flows are made available for the stream function and 
velocities due to distributions of circular vortices over axi -symmetric 


surfaces . 










Figure 3. Optimum Radial Distributions of Circ ula tion 

















Figure 8. Effect of Second Order "Swirl” Terms on the Streamtube Shape for 




Figure 9. Stream-tube Shapes for Various Values of the Thrust Coefficient and a Comparison 
of the Optinma and Piecewise Constant Axial Velocities in the Ultimate Wa>e for 





Figure 10. Streamtube Shapes for Optimum and Off -Optimum Circulation Distributions 
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APPENDIX A 


EXPLICIT FORMUIAE FOR AXI-SYMMETRIC VORTICES 

In studies of flow fields induced by rotors and propellers, 
information about velocity components and stream functions associated 
with various axi-symmetric distributions of vortices is frequently 
required. In the course of the present study, explicit formulae ■ p ->r 
the velocity components and stream functions have been developed for 
several important distributions of axi-symmetric vortices. The availa- 
bility of these explicit formulae facilitates the computation of the flow 
field associated with rotors and propellers. Selected formulae considered 
to be of general interest are presented in this Appendix. Of equal, 
perhaps even greater, utility is the approach developed for the 
derivation of these explicit formulae. This approach can be easily 
used to obtain explicit formulae for many additional distributions of 
axi-symmetric vortices. The approach is described in this Appendix. 

It is worth noting that analogous problems involving axi-symmetric 
distributions of singularities exist in maiqy other fields of study such 
as electrostatics, magnetostatics, and potential flows about axi- 
symmetric bodies. The approach described here Is easily adapted to 
these analogous problems. For example, as an offspring of the present 
study, explicit formulae for the velocity components due to sources 
of uniform strength distributed over a circular cylindrical segment 
were obtained (Ref. 32) and utilized in conjunction with the A. M. 0. 
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Smith method (Ref. 33) for computing potential flows associated with 
lift-fan inflow. 

It is well known that formulae for axi -symmetric flows frequently 
involve complete elliptic integrals. Complete elliptic integrals of 
the first, the second, and the third kinds are respectively defined as 


TT 


K(k) = f 


r 

o (l - k 2 sin 2 a ) 2 


(A-l) 


TT 


E(k) = J 2 (1 - k 2 sin 2 a ) 2 bx 


(A-2) 


and 


n (n\k) = f 


dor 


o (1 - n sin^ cr)(l - k 2 sin 2 a ) 2 


(A-3) 


The behavior of the complete elliptic integrals are well under- 
stood. Tabulated values as well as computational formulae for these 
integrals are available in standard handbooks of mathematical functions. 
Explicit formulae in terms of the complete elliptic integrals are there- 
fore convenient formulae to employ in numerical procedures. Many defied ' 
integrals involving sines and/or cosines of the variable of integration 
are expressible in terms of the complete elliptic integrals. Eor 
example, the integral 


rr 


2 2 

( k ) _ J2 sin g cos d tt 


(1 - k sin a) 


07 = 
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(aln g - sin „1 d„ 


o (l - k sin a) 


may be reexpressed, upon integrating by parts, as 



1,00 = if ^ 

1 k 2 J o (l - k 2 Sin 2 <*) 2 
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da 


o (l - k sin a) 2 k‘ 


2 f 2 , ,2 .2 1 

J-,TJ 0 ( lJc sin a) da J 


Thus one has 


J2 (sin |^sin g]_|g = -^ [(2 - k 2 ) K(k) - 2E(k)J (A-4) 


.4 


o (l - k sin a) 


As another example, consider the integral 


J 2 = 


^ cos B dB 
o (a - b cos (3)^ 2 


By letting (3 = 2(<y + tt/2) and recognizing that the resulting integrand 
is symmetric about a = 0, one obtains 


TT 


I it p (2 sin « • 1) to 

2 (a + ^ 2j o (l-k 2 sin 2 .)^ 


(A-5) 


where 


*=(^f 


Rearranging the integrand of Eq. (A-5) gives 


l 2 m tf /2 ,2 / TT \} 2 " 2k + k ^ J 2 ~ ^2 2 072 

(a+b)'k(l-k) o (1 - k sin o) ' 
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- (2 - k 2 ) J 2 (l - k 2 sin 2 a ) 2 


da 
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ft - * £ > f ] 


(1 - k sin a) 


Using Eqs. (A-l) , (A-2 ) , and (A-4), one then obtains 



cos 0 dp 
(a - ft cos p)^^ 



- 2K(k) + ( 



(A-7) 
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Circular Vortex Filament 


The velocity dq due to sin element of vortex filament ds is given 
by the Biot-SavsLrt Law: 



(r ' - r ) x ds 



(A-8) 


where h is the strength of the vortex filament, r' is the position 
vector of the vortex filament, and r is the position vector of the 
observation point where dq is evaluated. 

In a cylindrical coordinate system (r, 0, z), with a symmetric 
distribution of vorticity about the z-axis, it is easy to show that the 
r and z velocity components, u and w, are dependent only on the 0- 
component of vorticity, T|. Thus the basic building block for the ccm- 
putation of u and w are the velocities and the stream functions associated 
with circular vortex filaments . These quantities associated with a 
circular vortex filaments shall be designated by the subscript "f". 

Consider a circular vortex filament of strength k and radius R 
centered about the z-axis in the z = Z plane. Since u^, and axe 
independent of 0, one obtains from Eq. (A-8) 


u f (r, z; R, Z) = HR ^-"- Z ^ 




cos 0’ d0' 


o [r^ + R^ - 2rR cos 0' + (z 


ww 


(A-9) 


and 


w f (r, z; a, Z) = jg J 


r 


(R - r cos 0' ) d0* 


[r*" + R^ - 2rR cos 0' + (z - Z)^] 


7W 


(A-10) 


Equation (A-9) yields, with the use of Eq. (A-7), 


u f (r, z; R, Z) = Jn r~6 ^ ^ “ K ^] (A" 11 ) 

where 

8 = [(r + R) 2 + (z - Z) 2 ]* (A-12) 

and i 

k = (A-13) 

By writing the integral in Eq. (A-10) in terms of the integrals 
in Eqs . (A-l), (A-2), and (A-4), one similarly obtains from Eq. (A-10) 

W (r, z; Rj Z) = [ 2r - - R -^ E(k) + K(k)1 (A-l4) 

f 2 tt 6 L 2r (x . k 2 ) J 

The stream function Y^ due to the circular vortex filament is 
given by 

Y f (r, z; R, Z) = J r'w f (r\ z; R, z) dr’ (A-15) 

where we have set Y^(o, z; R, Z) to be zero. 

Instead of placing Eq. (A-l4) into Eq. (A-15) and performing the 
integration with respect to r', it is simpler to use Eq. (A-10) and 
integrate first respect to r ' . Thus , one writes 
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Integrating the first double integral by parts with respect to 9' gives 
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One then obtains 
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Letting 9* = 2(<> + tt/2), one has 
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Therefore 
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* f (r, z; R, Z) s jjJ {(2 - k 2 ) K(k) - 2 E(k) } 


(A-19) 


The velocity components u and w associated with various axi-symmetri 
cally distributed vortices are expressible as integrals involving u f and w f 
Consider, as an example, a distribution of circular vortices of constant 
strength X over the surface of a circular cylinder of radius R centered 
about the z-axis and extending between the planes z = z^ and z = Zg. 

The vortices are in the 0-direction and the velocity components u and 
w due to this distribution of vortices are expressible as: 


and 


u(r, z) 


r 2 

X J u f (r, z; R, Z) dZ 
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w(r. 


z) = 



z; R, Z) dZ 


(A-20) 


(A -21) 


Attempts have been made in the past to place Eqs. (A-ll) and 
(A-l4) into integrals similar to those appearing in Eqs. (A-20) and 
(A-2l) and to obtain explicit formulae for u and w for various axi- 
symmetrically distributed vortices. Such attempts, however, have met 
only with limited success. In particular, very few formulae are avail- 
able in the literature for the velocity components and stream functions 
for general field points. In the present work, an alternative approach 
was utilized and many new expliciv formula were derived. 


J 


69 



Cylindrical Vortices 


Consider the uniform distribution of vortices over a finite 
circular cylindrical surface described earlier. The velocity components 
due to this distribution of vortices are, according to Eqs. (A- 9 ), (A-10), 
(A-20), and (A-21): 


u (r, z) = j (z - Z) dZ ^ 
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(A-23) 


Instead of using Eq. (A-ll) to express the right side of Eq. 
r 4.-22) in terms of an integral involving complete elliptic integrals, 
the integration for (A-22) is performed first with respect to Z, yielding, 


u(r > z) ■ sl[ J 
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(A-24) 


One therefore has, using Eq. (A-l8), 


u(r, 
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where 
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Similarly, integrating with respect to Z, one obtains from Eq. 


(A-23) 
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Noting that 
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one obtains 
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(A-31) 


Equations (A-21) and (A-27) are convenient explicit formulae to 
use in the computation of velocities induced by cylindrical vortex 
segments. In computing rotor induced flowfield, these formulae offer 
the possibility of representing the trailing vortex system in the slip- 
stream by a set of vortex cylinder segments. The strength rnd radius 
of the vortex segments may be determined in such a way that the kinematic 
(Y = constant along each vortex tube) and the kinetic (force free) con- 
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ditions are satisfied. 


The behavior of u and w as given by Eqs. (A-25) and (A-3l) can 
be determined by examining the properties of the complete elliptic 
integrals. In particular, it can be shewn that u is well behaved 
everywhere except at the ends of the cylinder, i.e. at r = R and z = 
z^ or z 0 , where it possesses a logarithmic singularity, w is well 
behaved except on the cylinder where it changes discontinuously by the 
amount X across the cylinder. 

Equations (A-25) and (A-3l) yield directly a number of velocity 
formulae for specialized cases. Some of these specialized formulae 
have been studied by other investigators (Ref. 26, 31). These specialized 
formulae include that of the velocities on the cylinder of radius R 
itself, the velocities on and away from the cylinder due to a semi- 
infinite cylinder (z^ = 0, Zg = ®), and that due to an infinite cylinder. 
The formulae presented here, i.e., Eqs. (A-25) and (A-31), can further 
be utilized to derive an explicit formula for the stream function. 

This formula for the stream function is also useful in axi-symmetric 
flowfield computations. 

Consider, for the moment, a semi-infinite vortex cylinder of 
strength X and radius R, centered about the z-axis and extending from 
z « z 1 to infinity. The ve3 •'city components induced by this semi- 
infinite vortex cylinder shall be designated u^ and w^. Equation 
(A-25) gives, with Zg = ®, kg = 0. Equation (A-25) therefore reduces 
to 
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n i R + r 


-i ,2 2 1 R - r 

o 1 - n sin a 


(A-33) 


Therefore, for r < R, Eq. (A-3l) gives 


w. (r, z) « 5 


XU - *!> ^ r r 


Utt ( rR ) * 


[ HI 11 ("Vjl) - *<V ] 


and, for z = z^ and r < R, one has 


w i(r, 2 i) " | 


(A-3*0 


The stream function for this semi-infinite vortex cylinder 
can now be evaluated by using the relation 


^(r, z)*j r’w(r', z 1 ) dr' - r J* u(r, z') dz' 


(A-3 5 ) 


where the va.’uf <-,f is taken to be zero along the z-axis. 

Placing Eqs. (A-32) and (A-3^) into Eq. (A-35) and evaluating 
the first integral gives, 


V r -*> - t + ^ f [ ( h • \ ) K <V> - S7 E( V>] dz ’ (A ‘ 36) 
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where 


Y-[ 


4rR 

(r + R) ? + ( 2 ’ - 


Z l> 


2 J 


Equation (A-37) gives 


(A-37) 


4rR dk. 


dz' 


k^' 2 [4rR - (r + R ) 2 


(A-38) 


Using Eqs . (A-l), (A-2), and (A- 38 ), Eq. (A- 36 ) is reexpressed as 

V r ' ■> - f (1 - 1 afo. 

' ' *n xfn - x )* o (1 - x sin a )* 


Integrating by parts with respect to a gives 


Y (r z ) . Xr! . X(rR) 3/ ^ ^ to [2 sin 2 ft .oo6 2 g 

Y l (r ’ } ^RTJ n o (1 - x sin 2 .)3> 


(A-39) 


Integrating with respect to x then yields 

y r r -) „ Xr 2 X(rR) ^ ^jl 7 - ~ 2 lt if fsU 2 a - sln 4 tt] dt, 

Yl ’ ‘ ^ * 2 ” J o (1 - n teT Y* .lA>* 

The integral in the above equation can be expressed in terms of the 
complete elliptic integrals. One obtains, af „er some algebraic mani- 
pulations , 

. 2 X(rR)^ k. |z - z. | r . . , 

*> ■ H- ♦ 5F — - [ ( -V ‘ir ) :; <Y 
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- A E <*i> + H 1 n ("V j 


(k-ho) 


The stream function Yg(r, z) due to a semi-infinite vortex 
cylinder of strength X and radius R extending from z = z g to infinity 
and centered about the z-axis is, from Eq. (A-4 0), obviously 

*2 


~r ECkg) + ~~ n (nVg) 1 

*2 


(A-41) 


By the use of the principle of superposition, one then obtains 
the following formula for the stream function ¥(r, z) due to a finite 
vortex cylinder of strength \ and radius R centered abet the z-r 
and extending between the planes z^ and z 


or 


T(r, z) = Y 1 (r, z) - Y 2 (r, z) 


T(r > z) = i &F{ k i> - 2 xl [ ( j 5 + K(k i ) 


' 4 E < k l> + H 1 " ,j 

k l 

*2^ - *2 1 [ + ]} (A-te) 

*2 *2 
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For cylindrical vortices with a non-uniform distribution of 
strength, explicit formulae for the velocity components and for the 
stream function can be obtained by utilizing the approach and the 
results presented above. For example, with a linear- distribution of 
vorticity strength, X, i.e., 


X = Xi + X 2 z 


(A-43) 


where X 1 and Xg are constants, the velocity components and the stream 
function are each expressible as a sum of two parts, one part due to 
X x and the other part due to XgZ. The part due to X-^ has already been 
expressed in terms of the complete elliptic integrals [Eqs . (A-25), 
(A-31), and (A-42)]. To demonstrate the procedure for establishing 
expressions for the part due to AgZ, consider the radial velocity u^. 


“i (r ’ 2) = inr. 


\ 2 R p Z 2 


Z(z - Z) dZ 


r 


cos 9' d9* 


[r 2 + R 2 - 2rRcos 6' + (z - Z) c } 




(A-44) 


Equation (A-44) may be rewritten as 


/ x X 2 R r 2 r 2 " (r 2 + R 2 - 2r R cos 0') cos 6* d0’ 

isr J «J 0 + . 2r R oos ; lE . ww 


X 2 Rz S2 r 

+ ~W J (z ' z) 32 . 


2tt 


cos 0' d0' 


0 [r c + R 2 - 2rR cos 0’ + (z - Z) 2 ] 3/2 


XqR r Z p p2TT 

1ST J 42 J .2 . „2 


cos 0* d0' 


2T 


z, o [r + R - 2rR cos 0' + (z - Z) ] 


(A-45) 
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The double integrals in Eq. (A-45) can be readily expressed in terms 
of the complete elliptic integrals using the approach presented earlier. 
In particular, the first double integral is easily expressed in terms 
of the complete elliptic integrals by first performing the integration 
with reject to Z. Tho second double integral is identical to that 
appearing in Eq. (A-22), which has been expressed in terms of the com- 
plete elliptic integrals [Eq. (A-25)]. The third double integral, is 
easily shown to be equivalent to the second integral appearing in Eq. 
(A-35), which again has been expressed in terms of the complete elliptic 
integrals . 


Annular Vortices 


Consider a distribution of circular vortices of strength X^/ *4 
where is a constant, over an annular surface in the z ■ Z plane. 

The surface extends from r = r^ to r = r^. 

The velocity components due to this annular vortex surface are, 
according to Eqs. (A-9) and (A-10), given by 


X 3 (z - Z) p r 2 ,2 tt 


cos 9* d8' 


r \ J ' U 

U(r ’ Z) = ““ 1 “ 4, J o [r 2 + H 2 - 2rK cos 77 ~ Z TF* 


(A- 46 ) 


X~ - 2 ^TT 
w(r, z) = ^ j dR J 

r l 


(R - r cos 9* ) d9* 

o [r c + R c - 2rR cos 8' + (z - Z) 2 ] 3 ^ 2 


(A- 47 ) 


Integrating with respect to R gives 


X-,(z-Z) p^TT 

u(r > i) ■ l.j 


cos 8’ (rg - r cos 8*) d8’ 


[r 2 +(z-Z) 2 - r 2 cos 2 0P] [r 2 + r^ 2 - 2rr ? cos(5+(z~Z) 2 ]^ 


r 


cos 8'(r^ - r cos 8') d8' 


o [r 2 + (z-Z) 2 - r 2 cos 2 0F] [r 2 + r^ 2 - 2 rr^ cos 8' + (z - Z) 2 ] 2 


T ] (A - 48 ) 

13 J 


X-. ,2tt 

w ( r > z > * - nJ 


J 


48 * 


o p p 1 — *' 9 ) 

o [r + R^ - 2 r R cos 8' + (z - Z) ] 2 


It is ea'y to show, by letting 0' ® 2(<* + tt/2) and using partial 
fractions, that 
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i -2 ™ cos Q 1 (R - r cos 6') d9' 

^ o [r 2 + (z - Z) 2 - r 2 cos 2 9'] [r 2 + R 2 - 2r R cos 0' + (z - Z) 2 ] £ 


j? (2 sin 2 cr - 1) (R + r-2rsin 2 ty)d tt 

\ 2 -, p - J _ /, _ , _ _• -.2 \ /, _ v 2 . 2 


(z-Z) 2 [ (r + R ) 2 + (z-Zf] a o (l-i&jSin a)(l + fflgSin cy)(l-k sin a ) z 


2 2 


Z) 2 [(r + R) 2 + (z 


Z ) 2 ] 2 


J 1L Z-ZL * 


r J 


2 _ . 2 


o (l - k sin a) 2 


r ly.-ZL - 1 p Jn 5 x 

®2 ^ o (1 - sin a)(l - k sin or ) 2 


fz - Z) c + R + r I r2 da 1 

2r m l -* o (1 + mg sin 2 <y)(T- - k 2 sin 2 a)‘ 


^ [r 2 + (z - Z) 2 ] 2 + r 


^ [r 2 + (z - Z) 2 ] 2 - r 


Using Eqs. (A-l) , (A-3), (A-44) and (A-48), one then obtains 
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\^(z - Z) - Tp p / 1*2 r 


m i ) n ^ m i\ k i 


+ r 


- ( 2 -—-**) " ( -"^ k r > 


“I 

j 


r, + r 


- — | [‘ ,K(k r, ) ' ( 2 ' ~^r "l ) " (n J> k r> 
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i* r 
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where 


4rr. 


1 

2 


k = | 
r, L 


1 L (r + r x ) 2 + (z - Z) 2 


and 


r 2 (r + r 0 T + (z - Z) 


4rr 2 _ 

2 J 


It is easy to show that Eq. (A-45) gives 


\~ k r, 


w(r, z) = ^ r -| K(k ) t K(k ' j 

2 nr 2 r * r 2 r, f T 1 J 
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(A-53) 


(a-54) 


(A- 55) 


(A-56) 
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According to Eq. (A- 17), the stream function due to the annular 
vortex surface is given by 


*( r » z ) = j 


\_,r „ 2 2-v. 


cos 0' d0' 


r ± * c [r 2 + R 2 - 2rR cos 9' + (z - Z) 2 ] ? 


(A-57) 


Integrating by parts with respect to 0' yields 


, \ V r 2 r 217 R sin 2 0' d0* 

y(r, z) = ! dR | — p 2 “ ) x 2 i/2 

4tt -x J o [r + R 2 - 2rR cos 9' + (z - Z) 2 ]^' 2 


A 0 r 2 : 2 t \ 


” %Tj 


'cos 0' rR - r 2 - (z-Z) 2 ~l sin 2 9' d9' 


[r 2 +(z-Z) 2 - r 2 cos 2 0'][r 2 + R 2 - 2rRcos0' + (z-Z) 2 ] 2 


(A-58) 


R = r. 


The last integral in Eq.‘ (A-58) can be re -expressed as 


16 p2 [2rR sin 2 g - rR - r 2 - (z - Z) 2 l(sin 2 a - sin* g) dc 

fi(z - Z) 2 ' o (l - k 2 sin a) 2 (l - m^sin a ) (l + sin a) 


{[R 2 - r 2 • (z - Z) 2 ] K(k) [(R + r) 2 + (z - Z) 2 ] E(k) 


( ^ + 1 - f ) n (m i Nk) + f ( q ' A + ? ) n( ^ k) 1 


Therefore, one obtains fiom Eq. (A-5S) 
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) n (-»-\k ri ) } ) 


(A-59) 


It is clear that the approach described here permits explicit 
formulae for the velocity components and stream function to be obtained 
for other types of distribution of vortices in annular regions. 
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APPENDIX B 


OPTIMUM PERFORMANCE OF ROTORS AND 


PROPELLERS IN AXIAL FLIGHT 


Equations relating the thrust and the power of rotors or pro- 
pellers in axial flight to the velocity com* • nents in the ultimate 
wake can be established rigorously by a simple extension of the analysis 
of Chapter IV. They are: 


f v w 


P = cO i v w r ds 
^ WWW w 


T = if j t(w w - « m f + (20 r w - v w ) T„] as w 


(B-2) 


where w ffl is the ascent speed of the rotor or the advance speed of the 
propeller. 

The relation between the axial and tangential velocity components 
in the ultimate wake, which also can be established rigorously by an 
extension of the analysis of Chapter IV, is: 


2 2 2 R 2 

w w v „ w v 

Hr - n T w r „ + Hr - Hr + J Hr ^ 

' 


(B-3) 


Following the approach described in Chapter IV, a perturbation 
Av w in the tangential velocity is introduced at r^ = a and the resulting 


changes in the axial velocity distribution, the power, and the thrust, 
are sought. It is obvious from Eq. (B-3) that the relation between 
the tangential velocity perturbation and the change in the axial 
velocity distribution is not altered by the presence of a non-zero 
speed of advance. Consequently, since Eq. (B-l) also is independent 
of w^, the expression for the change in power previously developed for 
the hovering rotor is applicable to the present case of a non-zero 
speed of advance. One therefore has 


AP = 2np0eAr w a 3 { (fla - v ya )+ 

wa 


wa 


~F 


2 , 

»a v r dr 

w w w | 

J o *w ; 


(BJO 


The change in thrust is dependent on w and is 


AT 


2 r 

= 2TrpeAr a l 2Qa - v 

W L i 


- (oa - v ) 
wa w _ wa' 

wa 


2v w 

wa » 


fei 


(B-5) 


a o w 


Equations (B-4) and (B-5) may be combined to express AP in the 

form 

AP = h(a) AT (B-6 ) 

where h(a) is a function only of a, the location where the tangential 


ti 
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velocity perturbation is introduced. Following the analysis of Chapter 
IV, the criterion for optimum performance is that h(a) is a constant, 
independent of a. In other words, the criterion for optimum performance 




P r w (rv + Nw ) rdr 

W oe 


27- - v„) 
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- r w - 
w w 


r v + Nw 
V W a 

w 

w 


w 


(n r K - v w )] 


(B-7) 


where N is a constant. 

It is easy to show that, upon setting w = 0, Eq. (B-7) reduces 
to the optimum performance criterion previously obtained for the hover- 
ing rotor. 

A differential form of the optimum performance criterion is 
obtained by di r erentiating both sides of Eq. (B-7) with respect to r^. 
The resu 1 .±ng equation relates the derivatives dv^/dr^ and dw ^/dr^. 
Diffe rentiating Eq. (B-3) with respect to r^ yields a second equation 
relating dv^/dr . and dw ^/dr^. After some algebraic manipulations, one 
then obtains the following two first order differential equations: 


dw 


— = 2w v (r - v )[2v t r - N(4r - v ) w 
~ www w /L w w ' v w' w 


~ 2 ~ 


dr 
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+ 2 ASr + 2v r 2 ]/F 
w W W J ' 
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dr 


— — s s v (3w^ 7 - 2 N*w ^ (3? - v ) 
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where 


or ^ ^ 3~ 2 r 2 rw /** \ 2 

F = w r - tii\w r + fANr - v r (r - 2v )Tw 
w w wv: L w ww'w v ' ^ V 


/ A^ AW A ^A **A \ld AV 

+ (An - v r ) v (r - v ) r 

v w w ' w w w' w 


(B-10) 


A = is a non-dimensional advance speed, v euid v* w are velocity 

* ~ 

components non-diir.ensionalized with reference to Or^, r w and N are non- 


dimensional and N with reference to 


At the rim of the ultimate wake, Eq. (B-3) gives 


W A 2 

O A . ~ / 0 'V \ 

—x — = — - + v (2 - v ) 
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(B-ll) 


Taking the upper limit of the integral in Eq. (B— 7 ) to be R w , 


one obtains 
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(B-12) 


Equations (B-8) through (B-12) permit the computation of optimum 
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distributions v 
w 

tivcly v and A, 
w 

o 

described in Ref, 
presented in Ref, 


and w w for any given values of N and A, 

An iterative method for computing v^ 

27 for various values of v . Sample 

o 

30. 
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or, alterna- 

and w. is 
w 

results are 
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